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Abstract 

We investigate the question of how the knowledge of sufficiently many 
local conservation laws for a model can be utilized to solve the model. We 
show that for models where the conservation laws can be written in one- 
sided forms, like dQs = 0, the problem can always be reduced to solving 
a closed system of ordinary differential equations. We investigate the Ai, 
A2, and B2 Toda field theories in considerable detail from this viewpoint. 
One of our findings is that there is in each case a transformation group 
intrinsic to the model. This group is built on a specific real form of the 
Lie algebra used to label the Toda field theory. It is the group of field 
transformations which leaves the conserved densities invariant. 



1 Introduction 



For a Hamiltonian system to be "exactly solvable" it must have a sufficient num- 
ber of conservation laws. Liouville proved that if a system with N degrees of 
freedom (i.e., with a 2A^-dimensional phase space) has independent conserved 
quantities with mutually vanishing Poisson brackets, then the system is inte- 
grable by quadratures [Q. This means that there exist (at least in principle) an 
algorithm for determining the state of system at time t from the initial data at 
time to = (say). 

There are (in two space-time dimensions) many non-linear field theories which 
are known to have infinitely many conserved quantities. Since there are many 
ways to count to infinity this fact alone is insufficient to conclude that the system 
is integrable; it is however a necessary condition. And, even if sufficiently many 
conservation laws are known, it is in general quite unclear how they can be used 
to construct the solution. The situation is rather the other way around, the 
existence and form of the conservation laws are deduced from an already known 
method of constructing the solution, like e.g. the inverse scattering method. 

In this paper we demonstrate how knowledge of conservation laws can be 
explicitly used to solve e.g. the initial value problem, for a particular class of field 
theories. 

1.1 Illustrative example 

The basic idea behind the method is well illustrated by the (utterly trival) ex- 
ample 

ddif = 0, (1) 

where d = d / dt + d / dx = dz and d = d / dt — d / dx = dz- Eq. @j has infinitely 
many conservation laws. However, they can all be generated from the two local 
laws, 

aj = 0, 9J=0, (2) 

where J = dip = 93' + tt and J = dip = 99' — tt, as the space integral of arbitrary 
polynomials V in J, dJ, d^J,..., or arbitrary polynomials V in J, dJ, d'^J,.... It 
follows from (H) and the ordinary rules of differentiation that 

5P( J, 9J, . . .) = 0, dP{J,dJ,d'^J,...) = 0. (3) 

Thus, all information about the infinite set of conserved quantities is already 
contained in J = J{x, t) and J = J{x, t). This matches the fact that to solve the 
problem we must determine two functions, ip{x,t) and 7i{x,t) = ip{x,t). 

At this point we should explain our notation. We do not consider a canon- 
ical formulation in light-cone coordinates, but rather at a fixed time t. All our 
conservation laws should thus be viewed as expressions in terms of the canonical 
fields ip and vr, with e.g. d'^ip being a short-hand notation for a more complicated 
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expression. The latter are found explicitly through repeated use of the equations 
of motion, i.e. in this example if = it, tt = ip". Here we find for n > 1 

d\ = {2d,r-' (vr + If') , d\ = {-2dS~' (vr - v') , (4) 

but with the non-linear equations to be considered later these substitutions be- 
come more complicated. 

Now, to exploit the conservation laws to solve the initial value problem, we 
first determine J(x,0) = Jo{x) and J(x, 0) = Jo{x) from the initial data ip{x,0), 
7r(x,0). Eq. (|]) implies a simple time evolution for J, J: 

J{x,t) = Jo{x + t), J{x,t) = Jo{x -t). (5) 

The fields if{x,t) and 7r(x, t) are now found by inverting their relations to the 
conservation laws, i.e. 

7r{x,t) = ^ Jq{x + t) + Jo{x - t) , Lp'{x,t) = ]^ Jo{x + t) - Jo{x - t) . (6) 

Since this is a first order ordinary differential equation for if the solution is de- 
termined only modulo an integration constant ipQ. This is an ambiguity which 
arises because the conserved currents J and J are invariant under the field trans- 
formation 

{ip{x, t), tt{x, t)) ^ {ip{x, t) + v^o, vr(a;, t)) . 

Thus, to find the complete solution we need an independent determination of say 
if{0,t). This can indeed be found by combining (^) with the equation of motion 
(|T]). Define q(t) = if>{0, t) and p(t) = 7r(0, t). They satisfy the equations of motion 

m=pit), pit) = ip"io,t). (7) 

In general this would not constitute a closed system of equations, since ip"{0,t) 
is a new unknown quantity. However, here it may be determined from the con- 
servation laws. It follows from that 

¥'"(0,t) = i[j^(t)- J^H)]. (8) 

The right hand side of this equation is a known, and may be used to close the 
system (0). 

The crucial points in the above analysis are (i) that we have a sufficient num- 
ber of conservation laws which can be written in "one-sided" forms like (H), so 
that their time development are easily found, and (ii) that the relations between 
the fields and the conserved densities are known and can be inverted. The first 
condition appears to be fulfilled for all (conformal) Toda field theories. In this 
paper we shall explicitly consider a few of the simplest cases. The latter relations 
turn out to form a set of (generally non-linear) ordinary differential equations. 
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whose order increases with the number of field components. Thus, the integra- 
tion of these equations introduces a set of undetermined integration constants, 
refiecting the fact that the set of conserved currents are invariant under a certain 
group of transformations on the fields of the model. To find the complete solution 
these integration constants must be determined from initial data. They play the 
role of angle variables in the theory. However, since there is only a (small) finite 
number of such undetermined quantities, most of the (infinite number of) angle 
variables of the model are also determined by the conserved densities. In this re- 
spect a one-sided local conservation law must contain infinitly more information 
than just the (globally) conserved quantities of the model. 

1.2 Contents and organization of paper 

The rest of this paper is organized as follows: In the next section we analyse the 
simplest [Ai) Toda field theory, i.e. the Liouville model, in considerable detail. 
We discuss in particular various methods for solving the initial value problem, but 
also show how the conservation of the energy-momentum tensor can be utilized in 
a systematic procedure for constructing the general solution. In fact, knowledge 
of the mechanism behind the general solution makes it easier to utilize for solving 
the initial value problem. We repeat this analysis in the following sections for two 
Toda field theories with two-component fields; in section 3 for the A2 model and 
in section 4 for the B2 model. Both of these models have one set of conservation 
laws in addition to conservation of the energy-momentum tensor. For A2 the 
additional tensor is of spin 3, and for B2 it is of spin 4. In both cases we show how 
the conservation laws can be used to derive general solutions of the equations of 
motion, and also indicate practical algorithms for utilizing these general solutions 
to solve the initial value problems. 

As indicated above, each Toda field theory is labeled by the name of a classical 
Lie algebra. Names were given originally by BogoyavlenskipJ while giving the Lax 
pairs for the models. This naming scheme is connected to the method first used by 
Leznov and Savelievp to construct general solutions to Toda field theories. Here 
the equations of motion are enforced through a zero curvature condition for a non- 
abelian gauge theory constructed from the Lie algebra in question. It is not clear 
from this formulation whether these Lie algebras have any intrinsic connection to 
the Toda equations. We show that a Lie group obtained by exponentiating one of 
the real forms of the complex Lie algebra is an intrinsic symmetry group for the 
system. The symmetry in question is the group of field transformations which 
leave the conserved densitites invariant. For the Ai Toda field theory this group 
is SL{2, R), for the A2 Toda field theory this group is SL{3, R), and for the B2 
Toda field theory this group is Sp{2,R). We conclude with a section describing 
the generalization to all Toda field theories. 

The Toda field theories have previously been considered by many authors. 
The treatment we feel is most similar to ours is by Bilal and GervaisQ, |5|. Apart 
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from additional detail, our presentation differs from theirs in viewpoint and fo- 
cus. Among other interesting discussions on Toda field theories we would like to 
mention the references[|, |[ 



2 The Liouville equation 

The Liouville model can be described by the Lagrangian density 

C=\{d^){d^)-ae-'^. (9) 

where a and a are arbitrary real parameters. One may in the classical case let 
a 1 and o" ^ ±1 after appropriate redefinitions of the fields. Only the case 
a = 1 gives a Hamiltonian which is unbounded from below. The equations of 
motion become 

ddip = -ae^. (10) 



2.1 Qualitative behaviour and averaged motion 

We expect the solutions of ([lOD to have very different behaviour when cr = 1 
and a = —1. Assume x-space to be a circle of unit circumference. Then (p!0|) 
can be thought to describe a string which is wrapped around a cylinder Si x R, 
parametrized by the coordinates {x,ip), and can slide along the cylinder. Thus, 
ip measures the position of the string along the cylinder axis. There is also an 
external potential cre'^ in this direction. To obtain a qualitative understanding 
of the behaviour of the system consider the average position 

= {^{■,t)) = f'dx^{x,t). 
Jo 

With the approximation {e'^^''^^) = e^'^^''^^\ which is exact when is independent 
of X, we get a closed equation of motion for $: 



This equation has the solutions 



9«2 

e*(*) = — (11) 

cosh^ [a{t — to)] 

when cr = 1, and 

e*W = — (12) 

sinh^ [a{t — to)] 

when a = —1. In the latter case the field $ always becomes singular at some 
finite time. 



4 



When (7 = 1 the general (physically interesting) solution can be written as a 
fluctuation correction to the average motion (|lTD 

ip{x,t) = ^{t) + 6ip{x,t). (13) 

Since $(t) — > — oo as |t| oo (which makes the exponential term in ([T0|) van- 
ish) it must be possible decompose Sip{x,t) into right- and left-moving waves at 
asymptotically early and asymptotically late times, 

6ip{x,t) ~ (p^{x + t) + i^^{x — t) as t ^ — cxD, 

6ip{x,t) ~ (/?+(x + t) + (^+(x — t) as t — > oo, (14) 

On a cylinder the functions ip± and (f± must be periodic. 



2.2 The solution of Liouville 

The solution to ([T0| ) was already found by Liouville 0. It can be written in the 
form 

.^(..t) _ 8F'ix + t) G'jx-t) ^^^^ 
F{x + t) + aG{x - t)] 

where F and G are arbitrary functions of their arguments. Since (^) involves two 
arbitrary functions it can be expected to be the most general solution. However, 
many choices of F or G lead to solutions that are unphysical. Also, it is not 
entirely trivial to find which F, G correspond to a given set of initial data, e.g. 
(p{x, 0) and 7r(x, 0) = 0(a;, 0). The x-independent solutions (pJ],p!^) correspond to 
the case when F{^) = exp[a(,^ — to)] and G(0 = exp[a(^ -|- to)]- The asymptotic 
decompositions (|l^) can easily be achieved by writing (cr = 1) 

F(a; + t) =e'^("+*)Fi(x + t), Gix - t) = e"^""-'^ Giix - t). (16) 



By inserting these expressions into (p^ ) we find that 

~ 8e^'''{aFi + F[){aGi + G[)/Gl as t ^ -oo, 
e^ - 8e-^^'{aFi + F[){aGi + G[)/F^ as t ^ oo. 

Thus, by choosing Fi and Gi as the periodic solutions to 

F[ + aFi = ae'^-, G[ + aGi = aGje"^- , (17) 

we find a simple mapping between the fields at asymptotically early and asymp- 
totically late times: 

= ¥^_ -2 logFi, ^+ = + 21ogGi. (18) 
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Since Fi resp. Gi are uniquely defined functionals of (p- resp. (p^, wliich can be 
written down explicitly by solving (p!7[) , eqs. (|T^) define a canonical transforma- 
tion which is the classical analog of the S-matrix. 

However, for more general given initial data at (say) time t = 0, it is not 
so clear how F and G should be determined in the most convenient way. One 



(admittedly ugly) solvable system is obtained by differentiating (|T5D with respect 
to t and X before setting t = 0, 

F" IF' 1 . G" 2cr G' 1 , , , , , 

with (H) itself clS cl further (consistent) condition 

F'G' = \e'^^ {F + aGf . (20) 

Note that F and G are not uniquely determined from these equations; one has 
the freedom of choosing 3 arbitrary initial conditions for (|1^,^). This is related 
to the fact that the representation (^) for p is invariant under a 3-parameter 
(Mobius) group of transformations on {F,G). This group is generated by the 
transformations 

(F, G)^{F + a,G + aa) , (F, G) ^ c (F, G) , (F, G) ^ {f~\ G-') , (21) 
and may be viewed as a gauge group for the system. 



2.3 Solution by direct use of conservation laws 

Liouville's solution ([l5|) represents a situation where the mapping from the fields 
{ip, Tc) to the separated solutions (F, G) is quite complicated, while the inverse 
mapping is simple and explicit. As indicated in the introduction, one may instead 
attempt a different route to the solution, utilizing the fact that the system has 
conserved currents with simple time development. In this case the mapping of 
fields to separated free wave solutions is simple and explicit, while the inverse 
mapping is complicated. We now investigate this second possibility. 

The basic conserved densities for (|T^) are the light-cone components of the 
conformally improved energy-momentum tensor, 

T = T,, = {dipf - 29V, T = T,, = (dp)^ - 2d\. (22) 

They satisfy dT = and dT = 0. And, as with the example of the introduction, 
any polynomials P(T, dT, d^T, . . .) or V{T, dT, d'^T, . . .) satisfy BV = dV = 0. 
The space integral of the corresponding densities will all commute with the Hamil- 
tonian, and thus be elements in the Lie algebra of the 'group of the Hamiltonian'. 
This group has sufficiently many mutually commuting elements to make the sys- 
tem integrable. However, all the important information is encoded in the two 
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fields T and T, so these are the quantities one should utilihze to find the solu- 
tion. 

Expressed by equal time canonical fields we find 
T = {Ti + ^'f -A{ti' + + 2a e^, T = (vr - i^'f + 4 (tt' - + 2a e'^. (23) 
There are several ways to read and utilize these equations: 

1. First they can be used to find e.g. T(a:;, 0) = Tq{x) and T(a;, 0) = Tq{x) 
from the initial data, ^q{x) = (p{x,0) and 7ro(a;) = 7r(a;, 0). From the 
conservation laws for T and T it follows that their time evolution is given 
as T{x, t) = To{x + t) and T{x, t) = To{x — t). 

2. With T{x, t) and T{x, t) known and t fixed, eq. (123|) define a set of ordinary 
differential equations in x for (p{x, t) and 7^{x, t). They are no more compli- 
cated than (p!9|j20| ). However, the complete solution to these equations re- 



quires three initial conditions to be found, say q{t) = (p{0, t), r(t) = (p'{0^ t), 
andp(t) = 7r(0,t). 



3. Eqs. ( P3|) may be viewed at fixed x (say x = 0) as a means of expressing 
(p"{0,t) and 7r'(0,t) in terms of {q,r,p), T, and T. This may be combined 
with (|ToD to find the equations of motion for {q,r,p). They become 

q = p, 

r = ^pr + l(T-T), (24) 
p = l{p^ + r^)-lae''-l{T + T) 

Here it is understood T, T are to be evaluated at a; = 0. 

In summary, to find the fields at some fixed later time, given the initial conditions, 
one has to solve two sets of ordinary differential equations in succession^, first (^) 
and then (|23|). This is still somewhat less efficient than starting from the known 
general solution ([Tsl). However, the advantage of this method is its systematic use 
of the known conservation laws. Thus, it has potential for generalization to more 
complicated equations where the general solution is unknown. This method will 
work as long as a sufficiently large set of local conservation laws can be exphcitly 
found. 



2.4 The general solution by "separation of variables" 

In the previous subsection we indicated how the initial value problem for the 
Liouville equation could be solved by systematic use of its known conservation 

"'^Or, more generally, one can solve one set of equations along an arbitrary shaped curve in 
space-time, starting from a point on the (spacelike) curve of initial data 
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laws. The procedure involved the solution of two sets of ordinary differential equa- 
tions in succession. In most cases these equations must be solved numerically. 
With present day technology this is in practice — and as a matter of principle — 
comparable to the fact that even most analytic solutions at some stage must be 
evaluated numerically on a computer. (In this respect we still consider the numer- 
ical solution of ordinary differential equations to be significantly different from 
the solution of partial differential equations, both because of the computational 
effort involved and the number of parameters in the general solution.) 

However, if we want to compute the fields for many different times, or derive 
asymptotic relations like ([T8|), there is still considerable advantage in working with 
the general solution (|I5|) . For this reason we here propose and test a method for 
constructing the general solution from the knowledge of the conservation laws. 
Our main motivation is of course not to rederive (p!5D , but to investigate an ansatz 
which can be generalized to the more complicated Toda field theories. 

It is preferable to work with the field U = e""^/^. The energy-momentum 
tensor ( P^ becomes 

T = 4^, T = 4^. (25) 

Thus, the dependency of U on z cancels in T, and its dependency on z cancels in 
T. This suggests that it should be possible to factor out the relevant dependencies 
on z and z explicitly. To this end we write 

U = U{z,z) = u{z) h{z, z) u{z) = uhu, (26) 

and make the ansatz that T is a functional of u only, and T is a functional of 
u only. This means that the expressions (|25|) for T resp. T should hold with 
the replacements U —>■ u resp. U ^ u. These are linear 2nd order ordinary 
differential equations for u and u. Inserting (|26|) into (p5D , and using the above 
ansatz imposes two relations on h: 

du'^dh = 0, du^dh = 0. (27) 

It is straightforward to integrate these equations. Define functions H = H{z) 
and H = H{z) such that dH = u~'^ and dH = u'"^. Then we obtain the 
representations 

h = oiQ + aiH, h = ao + aiH, (28) 

where ) and ). These two representations are consistent with 

each other if each is a linear combination of the functions 1 and H. Introduce 
the vectors 

H = il,H), H={l,H). (29) 

Then both the linear equations ( P?] ) are satisfied by the representation (the nor- 
malization is for later convenience) 

h = 2-1/2 (h-X-H') , (30) 
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where X is a 2 x 2 matrix of real constants. We remind the reader that our 
definitions of H and H imply that u = (dH)^^^'^ and u = (pH^ ■ We have 
thus found a complete solution ansatz for U = uhu, which now must be inserted 
into (|l^). This gives 

hddh-dhdh = adHdH, 
which after insertion of the representation (|30D gives the condition 

(H-X-H^) [h-X-h^) - [H-X-hf) (h-X-H^) = a, 

with vectors h = h = {0,1). On evaluating the left hand side this reduces to 

det X = a. (31) 

With this condition on X, the conservation laws together with the ansatz (^) 
have indeed lead us to the general solution of the Liouville equation: 

2dHdH , ^ 

e'^ = J-- • (32) 

Liouville's solution (^) correspond to the case that X = J ^ , F(x + 1) = 

E{z) , and G[x - t) = -E{z) , with z = {t + x)/2 and z= {t-x)/2. 

2.5 The transformation group SL{2, R) 

We have found a solution parametrized by two arbitrary functions in one variable, 
and (assume a = 1) a matrix X G SL{2, R). What is the role of X? Note that 
H resp. H is determined from u resp. u, hence indirectly via (p5D from T resp. 



T. Equations (psD are 2nd order equations for u or u, and 3rd order equations 
for H or H. Thus, the general solution involves 3 arbitrary integration constants 
for H, and 3 arbitrary integration constants for H. It is, for arbitrary T and 
T, not possible to find explicit expressions for u and u. However, assuming that 
some particular solutions Uq and Uq have been found, it follows by variation of 
parameters that the general solutions to (^) can be written as 

U = + b J Uq'^^ Uq, U = (^a + b J Uq. (33) 

These results applied to the definitions of H and H show that when some par- 
ticular solutions Hq and Hq have been found the general 3-parameter classes of 
solutions can be expressed explicitly as 

H, = ^U^^, H, = hlhlh^ (34) 

03 Hq + 04 63 Hq + 64 
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' ai 


Ci2 \ 




V as 





where we without loss of generahty may choose ai 04 — 02 0.3 = &i 64 — 62 &3 = 1- 
By inserting (0) into (|32D we find that the substitutions 

Hq Hi, Ho — * Hi, 

are equivalent to keeping H and H fixed, and instead transforming X as 

= B'XA. (35) 

Thus, there is an SL{2, R) x SL{2, R) group of transformations on the fields such 
that the energy- momentum tensor is kept invariant. This group is a product of 
an SL{2, R) group of "gauge transformations" which leaves the physical fields 
unchanged, X A = X, and an SL(2, R) group which transforms the physical 
fields without changing the energy-momentum tensor. For a fixed H and H the 
matrix X parametrizes the equivalence class of field configurations which gives 
rise to the same conserved densities T and T. 

The Liouville equation can be formulated as a zero curvature condition con- 
nected to the Lie algebra Ai (which has sl{2,R) as one of its real forms). The 
above discussion show that a real form of the same Lie algebra occur in connection 
with physical symmetries intrinsic to the Liouville equation. 



2.6 Revisiting the initial value problem. 

The solution (^2]) now leads to a clean way of solving the initial value problem, 
with say the initial data given at t = 0. 

1. From the initial data we compute the conserved densities T and T. These 
can be found as functions of x, which at t = is directly related to z = x/2 
and z = —x/2. 

2. We next determine the functions u and u by solving (^51) . This will most 
likely have to be done numerically. The equations are of 2nd order, thus we 
need two initial conditions for both u and u. It is convenient to pick one 
point (say x = 0) on the initial data curve, and require that 

u{0) = m(0) = 1, du{0) = du{0) = 0. (36) 

3. With u and u completely known we find H and H by integration. It is 
convenient to choose the integration constants such that 

H{0) = H{0) = 0. (37) 

It now follows from (|36|) and the definitions of H and H that 

dH{0) = dH{0) = 1. (38) 
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4. It now remains to determine the matrix X. By inserting 
we find0 



H) into 



d (e-^/2^ 
B (e-^/^j 



(39) 



which together with ( |3TD gives X. By computing dde 'fi/'^ ^ = 2 X 
we again verify that the equation of motion is consistent with ([31|). 



This provides a complete determination of all quantities in (|3^), which now may 
be used to compute the fields at arbitrary times. The connection between the ini- 
tial data and general solutions of the Liouville equation have also been considered 
by Papadopoulos and Spence[p!0| 



3 Solution of the A2 Toda field theory 

The A2 Toda field theory can be defined by the Lagrangian| 

C = B(j)i (901 + 902 502 + 902 90; 



(40) 



where the CTj's are parameters which may be set to ±1 after appropriate redef- 
initions of the fields. Only the case cxi = (T2 = 1 gives a Hamiltonian which is 
bounded from below. We shall assume this case in our further analysis. The 
equations of motion simplify if we introduce new fields = |(0i + 02 — 203) and 
X = I (201 — 02 — 03), in which case they become 



ddip 



ddx 



(41) 



The third field degree of freedom, $ = 01 + 02 + 03, decouples and satisfies 
(9(9$ = 0. We note that (HII) is symmetric under the interchange ip ^ x- The 
equations are consistent with the assumption = x = 0, in which case they 
reduce to the Liouville equation. 



3.1 General solution 

Equations (^iD allow two basic sets of conserved densities. With the notation 
that = d'^ip, Xn = 9"x, the set corresponding to the conformally improved 

^Note that the left hand sides of eqs. ( [39| ) can all be expressed in terms of the initial data 
ipo and ttq. 

^For a conventional normalization of the kinetic terms wc will from here on use the definitions 
d={d/dt + d/dx)/V^ = d, siiidd={d/dt-d/dx)/V2 = dg. 
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energy momentum tensor is 



T = - (^ixi + X1-V2- X2, (42) 

and the same expression for T with the replacement d ^ d. The other set is two 
spin 3 conserved densities [O, 



Q = 2v5iXi - 2v?iXi - 2v5iV52 + 2x1X2 + V1X2 - V2X1 + X3, (43) 

and the same expression for Q with the replacement d d. Note that (T, T) is 
even and {Q, Q) is odd under the ^ x symmetry. This is seen more explicitly 
if they are expressed by the even = (</? + x)/2 and odd ip = ip — x fields. With 
the notation that = d"'(j) and ipn = d"'ip, 

T = 0? - 202 + Q = (20? - 02)^1 - - 301^2 + i^3. (44) 

These conserved densities satisfy dT = dT = dQ = dQ = 0. Thus, their 
time evolution is explicitly known. It is convenient to express them by the new 
variables U = e''^ and V = e~^. With the notation that f/„ = d^U and Vn = d"'V 
they become 

f/2 1^2 UiV, 

n = ^3 U3 U1V2-U2V1 U2U1 V2V1 UiV? ulVi . . 

and the same equations for T and Q with the replacement d —>■ d. In solving 
the initial value problem these relations should be viewed as equations for the 
fields U and V, with (T,T) and {Q,Q) already known quantities. We want to 
write them in separated form. To this end we use the conserved density ( ^5] ) to 
eliminate V2 resp. U2 in (|46|). We find 



jj-T-^ = -{dT-Q), y-T^ = -{dT + Q), (47) 

which are seen to be linear equations for U and V. There is a similar set of 
equations for U amd V, obtained by making the replacement T —>■ T, Q —>■ Q, 
and a in (0). 

As for the Liouville model we try to factor out the relevant dependencies on 
z and z exphcitly. We write 

U = u{z)h{z,z)u{z), V = v{z)k{z,z)v{z), (48) 

and make the ansatz that T, Q are functionals of u and v only, and that T, Q 
are functionals of u and v only. Thus u, v and u, v are determined by ordinary 
differential equations; those for u, v by making the replacement {U, V) {u, v) 
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in (^7[), and the similar ones for -u, by making the replacement {U,V,T,Q,d) — > 
(m, V, T, Q, B) in (|47|) . We still have to determine h and /c from equations in both 
z and z, but these will be explicitly given, and turn out to be quite managable. 

By inserting the factorized expressions for the U and V into (0), and de- 
manding that T and Q should be the above mentioned functionals of u and v 
only, we obtain equations for h and k. They may be cast into the form 



d — d — dh = 0, d — d — dk = 0. (49) 

U V V u 

There is another set of equations for h and k, obtained by making the replacement 
{u, V, d) {u, V, d) in those above. It is straightforward to integrate (49). Define 
functions H = H{z) and K = K{z) such that dH = v/v? and OK = u/v^. Then 

h = ao + aiH + a2 J KdH, k = + hK + ^2 j H OK, (50) 

where all the coefficients may depend on z: = ao{z) etc. Likewise, to solve 
the ^-dependence of h and k we define functions H = H{z) and K = K{z) such 
that SB = v/v? and BR = v/v"^. Then 

h = ao + aiH + a2 J KBH, k = f3o + Pi K + (32 J H BK, (51) 

where all the coefficients may depend on z: ao = ao{z) etc. 

Equations (^,^) are consistent with each other if each of the coefficients a, 
is a linear combination of the functions 1, H, J K OH, and each of the coefficients 
Pi is a linear combination of the functions 1, K, J H OK. Define 

H={l,H,JKdH), K={l,K,jHdK), (52) 

H = {1,HJKBH), K = {1,K, J HBK). (53) 

Then we have found that the h and k must be of the form 

h = HXH\ k = KYK\ (54) 

where X and Y are 3x3 matrices of real constants. 

Now inserting the ansatz = — log{v h v) and x = ~ log(t) k v) into the Toda 
equations (0) gives the equations h Bdh — Bhdh = k BH OH and k Bdk — Bkdk = 
hBROK. To evaluate these expressions further we first compute 

dH = {0,1, K) dH = hdH, OK = {0,1,H) dK = kdK, (55) 

and the similar barred relations. It then follows from ( ^4|) that 

hBdh= (H-X-H*) (h-X-h^) (BHOh] 



dhdh = (h-X-H^) (H-X-h^) (OHdH 
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and the same with {h, H, H, H, h, H, h, X) {k, K, K, K ,k, K , k,Y). We 
obtain the algebraic equations^ 



(h X H 



hXh} 
k-Yk' 



h X 
kYK^ 



H X 

K Y k^ 



K Y 
H X 



(56) 
(57) 



Regarding the expressions as polynomials in the indeterminates if, J K dH 
(with / HdK = HK — J KdH), and the corresponding barred quantities, (p^) 
are algebraic equations for X and Y. Their solution is 



P fx*)"' P, 



det X = det 1" 



(58) 



where P is the matrix 

' 1 
of H, K and if, K imply that 



1 

-1 I . We remind the reader that our definitions 



u 



{dHf dK 



-1/3 



dH {dKf 



'1/3 



(59) 



and the corresponding barred relations. With these expressions, eq. (^J), and 
the restrictions (|58D, we have found the general solution to (^). This solution 
is expressed in terms of the 4 arbitrary functions, (if, K, H, K) and a S'L(3, R) 
matrix (say X, with Y then been determined as the corresponding representation 
matrix for X). 



3.2 The transformation group SL{3,R). 

The form of the solution we found above strongly indicates that there is a group 
structure behind it. As with the Liouville equation this can be interpreted as the 
group of field transformations which leaves the conservation laws invariant. It 
turns out that this group acts linearely on the fields h and k. First note that the 
conservation laws are related to the fields ip and x through two 3rd and two 2nd 
order differential equations; thus there will appear a large number of integration 
constants in the solution of these equations. For a proper counting we find it 
most safe to (mentally) express everything in terms of the canonical fields (p, ir^p, 
X, T^xi ^"^^ their spatial derivatives. The number of integrations constants then 
turns out to be equal to the sum of the orders (spins) of the conservation laws 
minus the number of fields, here 2x3 + 2x2 — 2 = 8. 

Consider now our definitions (^) of H and K. With u and v given there is 
a freedom to shift H and K by constants, and the integral / K dH also contains 
an undetermined integration constant (its lower limit of integration). However, 

''if one wants to solve the model (^0|) with tri, (T2 € {1,-1} one should at this point multiply 
the right hand side of ( p6| ) by 02 and the right hand side of ( |57| ) by ui. These modifications 
can be undone by the substitutions Y ^ and X uiX . 
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such changes are equivalent to keeping H and K fixed, and instead multiplying 
X and Y by appropriate matrices. To be specific, changing — > if + a is 
equivalent to multiplying X resp. Y from the right by 



e"^i"'= U 1 , e"^i"'= I 1 I. (60) 
\001/ \Oal/ 

Changing K ^ K -\- (3 is equivalent to multiplying X resp. Y from the right by 

(61) 

Changing / K dH J K dH + 7 is equivalent to multiplying X resp. Y by 

(62) 








There is a similar freedom to shift the corresponding barred quantities. These 
shifts are equivalent to keeping H and K unchanged, and instead multiplying X 
resp. Y from the left by the transpose of matrices corresponding to those above, 
i.e. exponentials of generators E^^^ which are the transpose of the generators f/'^ 
above. 

Our equations arc also (for given T and Q) invariant under scale transforma- 
tions: u ^ V ^ V, h{u,v} — >• e^ /i{e^M, e'^v}, k{u,v} — >• k{e^u,e^v}. 
Such transformations are equivalent to keeping H and K fixed, and multiplying 
X resp. Y from the right by 

e'^^^"'+^<' = I '0 e--'* y e''<'+'^<' = I e''-- I (63) 
\ e-'' / \ e-'' / 

There is a similar invariance under rescalings of u and v. This is equivalent to 
multiplications from the left by the corresponding diagonal matrices. The two 
sets of 8 generators i?] (i = 1,2, j = 1,...,3) are representation 

matrices for a basis of the complex Lie algebra A2. Upon exponentiation of these 
generators with real parameters, and further composition, they generate groups 
of real matrices, more specifically two copies of SL{3, R) which are (inequivalent) 
representations of each other. 

The transformation on integration constants we have considered above only 
corresponds to the multiplication of X and Y by lower triangular matrices from 
the right and upper triangular matrices from the left. This does not generate the 
full group structure, because there are additional integration constants arising 
when we determine {u,v,u,v) from the conservation laws. But a change u — >■ 
u' corresponds to a multiplication of h by the factor u/u'. Thus, as with the 
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Liouville model, we expect that the freedom of changing all integration constants 
are equivalent to multiplying X by arbitrary elements of SL{3, R) from both 
the left and the right. A SL{3, R) subgroup of these transformations are gauge 
transformations which leaves the physical fields unchanged. 

3.3 The initial value problem. 

Let us summarize how the initial value problem can be solved using the results 
above: 

1. From the initial data at time t = to one computes the conserved densities 
T, T, Q, Q. Note that this process requires a rewriting of the quantities in 
terms of the canonical fields and their spatial derivatives. This can be done 
with the use of the equations of motion, although the expressions become 
rather lengthy. 

2. From the knowledge of the conserved densities one determines {u,v,u,v}, 
using (|47|) with the appropriate replacements mentioned above. For most 
initial data these equations will have to be solved numerically. Note that 
we may choose a fixed set of initial conditions in this process. 

3. With further integrations one finds the vectors ^H,K,H,K^. We may 
also choose a fixed set of integration constants in this process. 

4. Finally X and the corresponding Y aare determined so that the solutions 
(p = — log{u h u) and x = ~ log{v k v) reproduce the initial values of the 
canonical fields (and a sufficient number of their spatial derivatives) at some 
initial point {tQ,xo). 

4 Solution of B2 Toda field theory 

The B2 Toda field theory can be defined by the Lagrangian 

£ = 901 901 + 902 902 - icTie'^'^^-'^^) - a2 e^'^' , (64) 

where the (jj's are parameters which may be set to ±1 after appropriate redefi- 
nitions. Only the case ui = (72 = 1 gives a Hamiltonian which is bounded from 
below. We shall assume this case in our further analysis. The equations of motion 
simplify if we introduce new fields (p = 20i and x = 0i + 02, in which case they 
become 

dd^p = -ae^^-^^, ddx = -ae-^+2>^. (65) 
These equations have no apparent symmetry. However, the ansatz 

</. = 20 + log3, x= ^0 + log(3/v^), (66) 
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is consistent with the equations, and leads to the Liouville equation for 0. 



4.1 General solution 

Equations (^) also have two sets of basic conserved densities. With the notation 
that = Xn = d"'x the set corresponding to the conformally improved 

energy momentum tensor is 

1 1 
T = -y:>l- ViXi + xl- ^V2- X2, (67) 

and the same expression for T with the replacement d d. Eq. (B^) admit no 



spin 3 conserved densities, but there is a set of spin 4 conserved densities|12 



Q = (^t-4(^?xi + 12<^?X?-16<^iX? + 8xt-2<^?<^2-4(^iXi<^2 

+3 - 8 X2 + 24 xi X2 - 16 x? X2 + 4 ^2 X2 + 4 xl (68) 
+2 v?i + 2 Xl V^s - 2 v^i X3 - 4 Xl X3 - + 2 X4, 

and the same expression for Q with the replacement d ^ d. These conserved 
densities satisfy dT = dT = dQ = dQ = 0. If we insert the ansatz (|66|) into these 
expressions the result must be expressible in terms of the conservation laws of 
the Liouville equation. Indeed we find T = | T^, and Q = ^ — f d'^T^, where 
Tl is the energy-momentum tensor ( ^2]) of the Liouville model. This provides a 
check of our expressions. 

It is convenient to introduce new variables U = and V = e~^. We obtain, 
with the notation that f/„ = d'^U and Vn = d'^V, 



1 Uo Vo f/i Vi 

2 U V UV ' ^ ' 

and a lengthy expression for Q. It is possible to use (^) to eliminate all explicit 
reference to V resp. U in the expression for Q. We get 

^-E^ + ]l^.^T^ + 2T^-2dT^ = -Q + d'T-AT' (70) 

f/f/2 2t/2 f/^f/2 U 2^ ^ ^ 

and 

Yl^2Ty^-2dT^ = T'--d'T--Q. (71) 

V V V 2 4^ ^ ^ 

The latter is seen to be a linear Ath order equation for V. 

As before we write U = u{z) h{z, z) u{z), V = v{z) k{z, z) v{z), and make 

the ansatz that T and Q are functionals of u and v only (and that T and Q 

are functionals of u and only). This means that u must satisfy (^) with the 

replacement U —>■ u, and f must satisfy (|7Tp with the replacement V —>■ v. With 
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these assumptions we find equations for h and k. Tlie equation for h can be 
written in the form 

(^d-d-d^dh) h + R = 0, (72) 
yv' u u j 

where R is rather comphcated. With the notation that hn = d"'h, Un = d^u, and 
Vn = d^v: 

1 / u'^ u \ 

R = —u^hih^ H — u^hl — uuihih2 + 2ul — 2uu2 f2 uiVi h\. (73) 

2 y V V J 

However, the derivative of R has a compact representation, 

dR=-(^d — d — d^dh]dh. (74) 
yv' u u J 

Thus, by differentiating (^) we find that h satisfies a 5th order hnear equation, 

^,2 „.2 2 

9-;r9 — 9 — 9^9/1 = 0. (75) 

The equation for k can be written as 

7i2 „.2 

9 — 9 — a — afc = 0. (76) 

u u 

It is now straigthforward to integrate ( |75D and (|76|). Define functions Hi and -f^i 
such that dHi = v"^ /v? and dKi = u/v"^. Then 

h = QiQ + ai Hi + a2 H2 + 03 -ffs + 04 -f^4, 

A; = Po + PiKi+P2K2 + P,K,, (77) 

where all the coefficients a^, /3j may depend on z: ao = do{z) etc. The functions 
Hi and Ki are given by the integrals 

H2 = j dHi JdKi= J Ki dHi, 

H3 = J dHi J dKi JdKi = ^ J KldHi, 

H, = JdHiJdKiJdKiJdHi = HiHs-^Hi, (78) 

K2 = J dKi JdHi= J Hi dKi = HiKi- H2, 

K3 = I dKiJdHiJ dKi= H2K1-2HS. 
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Here the various identities among the integrals are derived by (repeated) partial 
integrations, and depends on certain assumptions about the undetermined inte- 
gration constants (i.e., the lower integration limits). This is the maximum set 
of relations which may be derived between the integrals without making spe- 
cific assumptions about the functions u and v. It is also the minumum set 
required for the rest of the solution to go through. However, the identities 
still allow us to freely change the outer integration limits in the definitions of 
H2 and H^, and to shift Hi and Ki by independent constants. Note that we 
have the relation 2 H4 — 2 H^Hi + if| = 0. With the additional constraint that 
2 «4«o — 2 a^ai + = the expression (^) for h is also a solution of ([72|), which 
in fact was the equation we first solved to determine h. 

The functions h and k must also satisfy differential equations in z, the same 
as (|75|) and (|76D with the replacements (m, v, d) {u, v, B). This means that we 
also have the representation 

h = ao + ai Hi + a2 H2 + «3 + «4 H4, 

k = P0 + P1K1+P2K2 + P3KS, (79) 

where the Hi^s and Ki^s are given by the barred version of (^), and the coefficients 
ai, Pi may depend on z: ao = ao{z) etc. Equations (^,|7|) are consistent with 
each other if each of the coefficients ai is a linear combination of the functions 
1, Hi, H2, H3, H4, and each of the coefficients is a linear combination of the 
functions 1, Ki, K2, K^. Define (with coefficients choosen for later convenience) 

H = (1, Hi, 2V2 H2, 2 H,, 2H^), K= (1, Ki, K2, K^) (80) 

H = {l,Hi,2'/^H2,2Hs,2H^), K = {1, Ki, K2, K^). (81) 

Then all the linear equations are satisfied by the representations 

h = HXH\ k = KYK\ (82) 

where is a 5 x 5 matrix of real constants and is a 4 x 4 matrix of real 
constants. However, not all the representations ( ^2t) will satisfy the original Toda 
equations (|65|). To check our ansatz we first compute 

dH = (0, 1, 2^/2 Xi, Kl, Hi Kl -2H2K1 + 2 H3) dHi = hdHi, 

OK = {0,l,Hi,H2)dKi = hdKi, (83) 

and the corresponding barred quantities. Inserting ip = —\og{uhu) and x = 
— log{v k v) into (^) gives the equations h ddh—dh dh = k"^ dHi dHi and k ddk — 
dkdk = hdKi dKi. These are purely algebraic conditions^: 

H X H^) (h-X-h^) - (h-X-H^) (S-X-h^) = (K-Y-K^Y , (84) 



(X-r-K*) (fc-r-fc*) - (fc-r-K*) (K-r-fc*) = [H-X-H^) . (85) 



^If one wants to solve the model with cti, (T2 G {1,-1} one should at this point multiply 
the right hand side of (p3) by cji and the right hand side of (p5|) by 02- 
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Regarding the expressions as polynomials in the indeterminates Hi, H2, H^, Ki 
(and the corresponding barred quantities) these are sets of algebraic equations 
for the elements of X and Y. 



4.2 The transformation group Sp{2, R) 

The equations (p^j85D are nonlinear and involve many (5^ + 4^ = 41) variables. 
They are not easy to solve by brute force. However, with the ansatz that X and 
Y are diagonal a two-parameter class of solutions is straightforward to find: 



X 



diag (1 



e'^,e2'^-^, 



l,e^-2^e- 



diag(e'',i 



e^-^,e- 



(86) 



The two parameters have their origin in the symmetry that our equations for 



given T and Q are invariant under: u 



V, h{u, v} h{e^u, e^v}, 



k{u, v} — > e'' k{e^u, e'^v}. Such transformations are equivalent to keeping H and 
K fixed, and multiplying X resp. Y from the right by diagonal matrices. 



exp (^/i h[ 
exp 



V 



diag (1 
diag (e'^ 



e^,e2^-'^, 



Le'^-^^e" 



(87) 



There is a similar invariance under rescalings of u and v. This is equivalent to 
multiplications from the left by the corresponding diagonal matrices. 

From the solutions (pB]) we can easily find many more, since the lower in- 
tegration limits in the expressions ([78|) for Hi, H2, H3 and Ki were not fully 
specified. We are free to shift Hi Hi+a. This is equivalent to keeping H and 
K unchanged, and instead multiplying X resp. Y from the right by matrices 



(X) 



( 1 

a 



V 







1 

a 



\ 




1 / 



aF- 



( 



1 




V 




1 

a 







1 / 



We are free to shift Ki 



BF, 



(X) 



(Y) 



Ki + p. This is equivalent to keeping H and K 
unchanged, and instead multiplying X resp. Y from the right by matrices 

/I \ 

1 

(3V2 1 

(3y/2 1 

V 1 / 

We are free to shift if 2 —>■ H2 + j. This is equivalent to keeping H and K 
unchanged, and instead multiplying X resp. Y from the right by matrices 



/ 1 



V 






1 / 



(X) 



1 






V -1' 





1 


-7V2 



\ 






1 / 



■yF. 



(Y) 



( 





1 



7 






1 / 
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We are free to shift ^ — 6/2. This is equivalent to keeping H and K 
unchanged, and instead multiplying X resp. Y from the right by matrices 



(X) 



( 1 




-6 

\ 





1 
1 



-(5 1 / 





/ 1 











\ 







1 








, e 4 = 








1 

















1 





There is a similar freedom to shift the corresponding barred quantities. These 
shifts are equivalent to keeping H and K unchanged, and instead multiplying X 
resp. Y from the left by the transpose of matrices corresponding to those above. 



i.e. exponentials of generators e\ ^ which are the transpose of the generators 
above. The two sets of 10 generators ^^hW f'-'^^ = 1, 2, j = 1, . . . , 4) 
are representation matrices for a basis of the complex Lie algebra B2. Upon 
exponentiation of these generators with real parameters, and further composition, 
they generate groups of real matrices. More specifically, with a symplectic form 
e and a metric form 77, 



(■) 



/ 1 ^ 

0-10 

10 

\ -1 y 



/ 1 \ 

-1 

10 

0-10 

V 1 y 



(88) 



the 4-dimensional representation generate real symplectic matrices S G 5*^(2, R) 
{S^ e S = e), and the 5-dimensional representation generate pseudo-orthogonal 
matrices O G 0(3, 2) (O* r)0 = rf). It is now straightforward to verify by direct 
insertion that a solution to ( p4| , ^5D is obtained by choosing any Y e 5*^(2, R) 



(as defined above) and taking X to be the corresponding 50(3, 2) representation 
matrix. The precise connection between corresponding X and 1^'s is implicitly 
defined by the connection between the 4- and 5-dimensional generators listed in 
the appendix. A more direct relation between X and Y is defined by (|55|), which 
moreover is an explicit expression for the function h. Eq. (|85D show that Y and 
—Y leads to the same solution for X, which means that Sp{2,R) is a double 
covering of 5*0(3, 2). Thus the transformation group in this case is Sp{2, R) (or 
rather Sp(2, R) x Sp{2, R) if we include the gauge transformations). 

The initial value problem can now be solved by the same procedure as de- 
scribed for the A2 Toda field theory. The practical implementation of this pro- 
cedure only becomes more cumbersome. 



5 Concluding remarks 

In this paper we have in detail discussed the explicit solutions of the three simplest 
Toda field theories. From this a general picture emerges, which seems to be as 



21 



follows: For each Lie algebra An {n > 1), Bn (n > 2), C„ (n > 3), Dn (n > 3), 
En {n = 6,7,8), F4, and G2 there exists a Toda field theory. The number of 
independent field components in the model is equal to the index n. Thus, to 
solve the model by the method used in this paper we have to find 2n independent 
conservation laws which can be written in one-sided forms, 

BQs = 0, dQ, = 0. (89) 

Such conservation laws can be found in all the cases we have considered. The 
conserved densities always occur in parity related pairs {Qs,Qs)- For there 
is one independent density for each spin s = 2,3, ...,n + 1. For S„ and 
Cn there is one independent density Qg for each spin s = 2,4, .. . , 2n. For Dn 
there is one independent density Qg for each spin s = 2,4,...,2n — 2, and an 
additional one of spin s — n. For Eq there is one independent density for each 
spin s e {2, 5, 6, 8, 9, 12}. For E^ there is one independent density for each spin 
s G {2, 6, 8, 10, 12, 14, 18}. For £^8 there is one independent density for each spin 
s G {2, 8, 12, 14, 18, 20, 24, 30}. For F4 there is one independent density for each 
spin s G {2, 6, 8, 12}. For G2 there is one independent density for each spin 
s G {2, 6}. Counting shows that the number of independent conserved densities 
is always equal to the number of field components. 

The fields in the model are related to a conserved density by a differential 
relation of the same order as its spin. By considering the canonical degrees of 
freedom one finds that the number N of integration constants which occur when 
one attempts to determine the canonical fields from the conserved densities is 
equal to the sum of the spins (counting both Qg and Q^), minus the number of field 
components. This number is iV = (n + 1)^ — 1 for An, N = (2n + l)n for Bn and 
Cn, N ={2n- l)n for Dn, and = (78, 133, 248, 52, 14) for (Eg, Et, Es, F4, G2). 
These numbers are precisely equal to the dimensions of the Lie algebras which 
label the models. Thus there is a A^-dimensional manifold of field configurations 
which lead to the same conserved densities, and a continuous transformation 
group acting on this manifold. Our experience from the previous sections is that 
the manifold may be identified with the transformation group, and that the group 
has as its Lie algebra one of the real forms of the (complex) Lie algebra which 
labels the Toda field theory. The freedom of choosing integration constants is 
probably such that the transformation group will consist of one simply connected 
component. 

How useful are the solutions we have found? If the purpose is to determine 
the classical fields on a dense grid of points covering a finite regular region of 
space-time, it may in fact be difficult to beat a direct numerical solution of 
the partial differential equations in efficiency (although probably in accuracy). 
The computational effort must in any case be proportional to the number of 
grid points, and the proportionality factor is rather small for a direct numerical 
method. However, if we are interested in finding the fields along a line or on small 
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set of points only, the use of our solutions would be vastly more efficient than a 
direct numerical method. The same is true if the interest is in how asymptotic 
field configurations at t — — oo develop into asymptotic field configurations at 
t — > cxo. Finally there is the problem of understanding the quantized versions of 
these models. For this any improved understanding of their analytic structure is 
of potential use. 



A The complex Lie algebra B^. 

The complex Lie algebra B2 can be represented by 5 x 5 or 4 x 4 matrices. 
Let diag"^ (•,...,•) denote a matrix whose only non-zero entries are on the n'th 
superdiagonal, and diag"^ (-,...,■) denote a matrix whose only non-zero entries 
are on the n'th subdiagonal. Then the transformations found in section 4.2 
correspond to a 5-dimensional matrix representation of a basis of S2, 

Hf^ = diag (1, -1, 0, 1, -1) , Hf^ = diag (0, 2, 0, -2, 0) , 

fP = diag^T 0, 0, 1) , Ef^ = diag^^ (0, 72, ^2, o) , 

Ff = diag^T (v^, 0, -x/2) , Ff = diag^^ (-1, -1) , (90) 

f[''^ = diagi^ (1, 0, 0, 1) , Flf^ = diagi^ (o, v^, v^, o) , 
= diag^i (x/2, 0, -x/2) , = diag^^ (-1, -1) , 

and a corresponding 4-dimensional representation. 



(91) 



These generators arc given in the Chevalley basis. Thus they are normalized so 
that [Fj, Fj\ = 6ij Hi for i,j G {1, 2}. Not all relative signs are determined by the 
conventions of a Chevalley basis. Here we have chosen signs so that [Fi, F2] = F3 
and [F2,F3] = 2F4. 





= diag (0,1, -1,0), 




= diag (1,-1, 1,-1), 


Ff) 


= diag^T (0,1,0), 




= diag^T (1 0, 1) , 


Ff 


= diag2T(-l,l), 


Ff) 


= diag^'T (1) , 


^(4) 


= diag^^ (0,1,0), 


Ft' 


= diag^^ (1,0,1), 


Ft' 


= diag2-^ (-1,1), 


Ft' 


= diag^^ (1) . 
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